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ABSTRACT 

The high-pressure carbon monoxide (HiPco) process for producing single-wall 
carbon nanotubes (SWNT) uses iron pentacarbonyl as the source of iron for catalyzing 
the Boudouard reaction. Attempts using nickel tetracarbonyl led to no production of 
SWNTs. This paper discusses simulations at a constant condition of 1300 K and 30 atm 
in which the chemical rate equations are solved for different reaction schemes. A lumped 
cluster model is developed to limit the number of species in the models, yet it includes 
fairly large clusters. Reaction rate coefficients in these schemes are based on bond 
energies of iron and nickel species and on estimates of chemical rates for formation of 
SWNTs. SWNT growth is measured by the co-formation of CO 2 . It is shown that the 
production of CO 2 is significantly greater for FeCO due to its lower bond energy as 
compared with that of NiCO. It is also shown that the dissociation and evaporation rates 
of atoms from small metal clusters have a significant effect on CO 2 production. A high 
rate of evaporation leads to a smaller number of metal clusters available to catalyze the 
Boudouard reaction. This suggests that if CO reacts with metal clusters and removes 
atoms from them by forming MeCO, this has the effect of enhancing the evaporation rate 
and reducing SWNT production. The study also investigates some other reactions in the 


model that have a less dramatic influence. 
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INTRODUCTION 

In the high-pressure carbon monoxide (HiPco) process the Boudouard reaction, 
catalyzed by iron, produces single-wall carbon nanotubes (SWNT). 1 ’ 2 Whereas using iron 
carbonyl as a catalyst yields significant amounts of SWNTs, attempts to use nickel as a 
catalyst have not led to SWNT production. At first, differences in the nucleation rates of 
iron and nickel were thought to have a significant effect on production. However, 
parametric studies of nucleation did not show a significant influence on production since 
there are other pathways to produce dimers and larger clusters without direct Fe + Fe 
nucleation. It was noted that one difference between iron and nickel is their different 
affinities for bonding with CO. Sunderlin, et al. ' report the bond energy for NiCO is 
about 170 kJ/mol (20468 K) and that of FeCO of about 35 kJ/mol (4214 K). The higher 
binding energy of NiCO will result in slowing down of its dissociation and possibly 
speed up its formation compared with FeCO. In the models we have included cluster 
growth by the exchange reaction 

Me n + MeCO => Men +1 + CO (1) 

where the reaction rate is taken to be the gas kinetic rate (Me=metal). These reactions are 
less likely for nickel because the binding energy of NiCO (170 kJ/mol) is only slightly 
smaller than that of Ni n -iNi (E a = 203 kJ/mol, n=2). This is not the case for Fe n .iFe where 
FeCO bond energy E a is 35 kJ/mol, whereas that of Fe-Fe is about two times greater, 75 
kJ/mol. The reverse of these reactions is the “CO enhanced evaporation” reactions in 
which CO abstracts Me from the cluster. Since the binding energy of NiCO is much 
stronger than FeCO we might expect greater loss of nickel clusters compared with iron 
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clusters. Unfortunately, we do not have any reliable rate coefficients for these reactions 
for n>0. 

Another factor to be considered is the direct or normal evaporation of small metal 
clusters, such as 

Ni n => Ni + Ni n ., (2) 

If they evaporate too fast, clusters may not be able to form, hence reducing the 
possibility of SWNT production. In the present models the smallest metal cluster that 
produces SWNTs is that of Feio or Niio. Therefore, clusters must grow to at least that size 
before SWNTs can be catalyzed. Several sets of rate coefficients are investigated for 
these evaporation reactions. 

These possibilities were studied parametrically using estimates of reaction rate 
coefficients based on limited knowledge of bond energies from the literature. The 
following discussion addresses the source of rate coefficients for various reactions in the 
models. These reactions are categorized as metal carbonyl dissociation and 
recombination reactions, dimer dissociation, and cluster evaporation (2), dimer 
nucleation, to some extent the exchange reactions (1), and the reaction 

MeCO + MeCO => Me 2 + 2 CO (3) 

Other reactions in the model, such as cluster growth and CO attachment reactions have 
been discussed in Ref. 4. 

Cluster evaporation data for iron and nickel can be estimated using the formula of 
Rao, et al. 3 (denoted “Girshick” in this article). 

The expressions for cluster growth /C and evaporation Ej are, respectively. 
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where © = — - is the dimensionless surface energy, and Vi is the hard-sphere collision 

kT 

frequency of monomers, p p is the density of bulk iron, i and j are the number of atoms in 
the colliding clusters, cr is the surface tension of the bulk metal. Si is the surface area of 
the monomer (atom), k is the Boltzmann constant, and T is the temperature. The 
saturation number density n s is a function of temperature determined from a Clausius- 
Claperon type relation, p s = p vre f exp(-E v /kT), where pvref is a reference pressure and E v is 
the energy required to remove an atom from the surface of the bulk substance. In the 
present study the saturation vapor pressure was determined by curve fits to the data of 
Ref. 6. Table 1 gives the parameters that were used in (4) and (5). 


RESULTS 

There are two sets of comparisons that were considered in comparing the 
influence of various reactions on the predicted growth of nanotubes. One is for iron and 
the other for a “cooked-up” model of nickel. The basic rate coefficients for the model are 
those of the Ames 2001 version 7 which were based on Krestinin, et al. 8 and calculations 
of rate coefficients for Fe(CO) x decomposition from Sunderlin, et al. 1 bond energies. The 
basic differences in these iron and nickel models are the rate coefficients for carbonyl 
reactions and the metal cluster growth and evaporation reactions. The basic model is a 
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variation of that described in Ref. 7, in which the iron clusters Fe n and their 
corresponding Fe n CO and CNT n clusters are lumped into groups having 8, 16, 32, 64, 
2048 iron atoms. Thus that model accounts for the agglomeration of clusters. Variations 
in the reaction rates that were assessed are those mentioned in the previous section. 
Definitions of these reactions and rate coefficient sources are given in the following 
paragraph. 

Carbonyl rates are defined by those reactions of the form 

Me(CO) x O Me(CO) x _i + CO (6) 

For iron, these rates were determined by Ames Research Center using the bond energies 
of Sunderlin, et al. ' For nickel, the rates were determined by simple substitution of the 
bond energy into the Arrhenius expressions for iron. No attempt was made to adjust the 
pre-exponential factors. 

For both the iron and nickel models the reaction MeCO + MeCO => Me 2 + 2 CO 
was taken from Krestinin original set of reactions with its rate coefficient for iron. This 
reaction is usually negligible. It only has an influence on the results when there is no 
other mechanism to form Me 2 . 

The MeCO bond energy exhibits its greatest influence when considering its 
dissociation. Two estimates were compared, one taken for the value of iron and the other 
for the value of nickel, 23.9 kJ/mole (2875 K Ames) and 170 kJ/mole (20433 K) based on 
Sunderlin’s data, respectively. 

Dissociation of metal dimers has a small influence on the results. The Krestinin 8 
value of the rate coefficient was used. These results and those using an estimate based on 
Girshick’s evaporation expression (5) were compared for iron. The surface tension in that 
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expression was adjusted until the activation energy matched that of the bond energy for 
the iron dimer 75 kJ/mol (or 9047 K) 9 . For nickel, the rate coefficient for iron was used, 
except the dissociation energy for Ni 2 was substituted in the Arrhenius energy term. The 
other method used expression (5) with the values for nickel substituted, 203 kJ/mol (or 
24476 K). 

Three models of evaporation of iron atoms from clusters were used. The first was 
the original set from Krestinin, et al. for iron. The second was based on (5) using the rates 
determined from the surface energy of Ref. 4, and the third was based on adjusting the 
surface energy until the activation energy for the dimer equals the bond energy of 75 
kJ/mole (9047 K). For nickel, three variations were studied: Krestinin*' values 
corresponding to iron, the Ref. 4 Girshick values, and the Girshick value calculated from 
(5) for nickel property data of 203 kJ/mole (24476 K) bond energy and its vapor pressure. 
Fig. 1 gives the evaporation rate coefficients for each model, evaluated at 1300 K. There 
are three basic sources of rates: Krestinin, et al. 8 , Rao/Girshick 4 equation (5) evaluated 
using nominal properties of Fe, and equation (5) with the surface tension adjusted so that 
the activation energy in the Arrhenius expression for dimer dissociation equals the value 
given in Ref. 9. 

The conditions of this study are for a constant temperature of 1300 K and 30 
atmospheres pressure, arbitrarily run for 0.1 seconds. The amount (mole fraction) of 
carbon dioxide produced at the end of this time was a measure of the effectiveness of the 
model. The starting amount of iron pentacarbonyl is 17 ppm in carbon monoxide. The 
basic cluster model is that called the Binary 2048 model. Metal clusters have 1-8, 16, 32, 
64, 128, 256, 512, 1024, and 2048 atoms. In this model there is a smaller number of 
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species so that the cases run very fast. The chemical rate equations were solved using 
program SENKIN of the CHEMKIN 3.61 package. The results of a set of calculations are 
presented in Tables 2 and 3 for iron and nickel, respectively. 

Iron 

The nominal iron case is denoted AD in Table 2. All the rates are based on the 
Ames model except that evaporation of the large iron clusters was determined from (5), 
Girshick. This model produced a CO 2 mole fraction of 0.001 in 0.1 s. Variations on this 
model using other values of Fe 2 dissociation and Girshick evaporation rates for small Fe- 
clusters only slightly affected the results. Reducing the Fe-Fe nucleation rate to zero 
reduced production by almost two orders of magnitude for the nominal case, but slightly 
increased CO 2 in the other cases where Girshick evaporation rates were used along with a 
smaller bond energy of 75 kJ/mole. When the smaller bond energy is used along with 

Q 

Krestinin, evaporation rates for small clusters, the amount of CO 2 is significantly 
reduced. It appears that evaporation takes its toll on cluster formation, hence their ability 
to catalyze the Boudouard reaction. The situation is exacerbated when the evaporation of 
clusters is based on calculations of the rate coefficients using (5) and the 75-kJ/mole bond 
energy, especially when the dimer dissociation rate coefficient uses this bond energy in 
its Arrhenius coefficient. 

To test the influence of the exchange of Fe from reactions with FeCO they we 
deleted in the last case in Table 2. There was a negligible effect on the results, indicating 
that these reactions are not important. 
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Nickel 

The models for nickel as a catalyst are based on more approximate estimates of the rate 
coefficients than the basic iron model. Probably, the most representative nickel models 
are cases no. 3 and 6 of Table 3. In those cases the rate coefficients have been modified to 
account for the bond energies associated with nickel. At 0.1 seconds the amount of CO 2 
produced in those two cases is S^XIO' 1 " and 9.92X10' 25 , respectively. This is negligible 
production, and reflects the nature of the strong NiCO bond 1 70 kJ/mole (20422 K). This 
strong bond prevents the formation of Ni clusters, thus inhibiting SWNT growth. The 
influence of the bond energy is demonstrated in cases 1 , 2, 4, and 7, where the NiCO 
dissociation rate was made the same as for iron. In those cases the production of CO 2 is 
similar to that for the iron model. Behavior similar to that of iron is seen in the influence 
of Ni-Ni dimer formation. When NiCO bond energy of 170 kJ/mole is used there is some 
influence of the NiCO + NiCO reaction and the NiCO + Ni n exchange reactions, as seen 
in cases 6, 8, 9, and 10. These reactions are neglected in the models of Dateo, et al. 

(Ref. 5). 

CONCLUSIONS 

Simulations of reactions in the HiPco process for SWNT production were carried 
out to investigate the influence of various possible reactions in the scheme and the effect 
of bond energy differences between iron and nickel. Time dependent calculations were 
made for a constant temperature of 1300 K and pressure of 30 atmospheres using a 
premixed reactor solution to the chemical rate equations. The principle metric for 
evaluating the models was the mole fraction of CO 2 produced by the Boudouard reaction, 
which is a measure of the total production of SWNTs. 
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It was found that the metal-CO bond energy has a major influence on the 
production of CO 2 . The higher NiCO bond energy almost eliminates production of CO 2 , 
whereas, the smaller bond energy of FeCO leads to significant production. Also, the rate 
of dimer dissociation and small cluster evaporation affects the production of CO 2 by 
limiting the rate of cluster growth, and thus the number of clusters available to catalyze 
SWNT growth. This suggests that if CO reacts with metal clusters and removes atoms 
from them by forming MeCO, and then this has the effect of enhancing the evaporation 
rate and reducing SWNT production. 

Because of the possibility of producing dimers via FeCO reactions, the rate of 
direct Fe + Fe => Fe 2 dimer formation did not have a big influence on CO 2 production. 
The same was observed for nickel if the bond energy of Ni-CO was assumed to be equal 
to that of iron, a nonrealistic circumstance. 
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Table 1 Parameters used in Equation (5) for cluster evaporation rates 


Model 

Surface 
Tension, cr 

N/m 

Ref. Vapor 
Pressure, P, T cf 
Pa 

E v 

Vapor 

K 

Atom 
Radius, ri 
M 

Atom Surface 

Area, Si 
2 

m 

Girshick 

Iron 

1.7 

6.83X10 11 

47400 

1.40X10' 10 

2.46E-19 

Girshick 
Fit to Ni 
Dimer bond 

2.53* 

9.64X10 11 

49170 

1.35X10" 10 

2.29E-19 

Girshick 
fit to Fe 
dimer bond 

3.65* 

6.83X10" 

47400 

1.40X1 O' 10 

2.46E-19 


^Inferred by adjusting it in Eqn. (5) to make E a equal to dimer bond energy 
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Table 2 Production of SWNTs with iron catalyst as indicated by C0 2 mole fraction for ten combinations of rate coefficients 
and dimer nucleation 


Fe Cluster Model T=1300 K P=30 atm 

Desig. 

Carbonyl 

Rates 

FeCO + 
Fen 

Fe2 

Dissoc. 

Fen Evap. 
8>n>2 

Fen Euap. 
na3 

C02 

Produced 
at 0.1 s No 
Fe-Fe 
Nucleation 

C02 

Produced 
at 0.1 s 
Fe-Fe 
Nucleation 
4E15 

AGd 

Ames 

yes 

Krestinin 

Ea=134 

kj/mole 

Girshick 

Girshick 

1.05E-03 

4.84E-04 

AG 

Ames 

yes 

Girshick 

Girshick 

Girshick 

1.05E-03 

4.88E-04 

AGdw 

Ames 

yes 

Ames 

Ea=75.2 

kJ/mole 

Girshick 

Girshick 

1.87E-05 

1.02E-03 

AD 

Ames 

yes 

Krestinin 

Ea=134 

kj/mole 

Krestinin 

Girshick 

2.37E-05 

1.00E-03 

AB 

Ames 

yes 

Ames 

Ea=75.2 

kj/mole 

Krestinin 

Girshick 

8.97E-08 

5.33E-06 

AW 

Ames 

yes 

Krestinin 

Ea=134 

kj/mole 

Eqn. (5) 
w/Ea=75.2 
kJ/mole 
bond 

Eqn. (5) 
w/Ea=75.2 
kJ/mole 
bond 

1.56E-18 

9.19E-15 

AWd 

Ames 

yes 

Eqn. (5) 
w/Ea=75.2 
kj/mole 
bond 

Eqn. (5) 
w/Ea=75.2 
kJ/mole 
bond 

Eqn. (5) 
w/Ea=75.2 
kJ/mole 
bond 

9.19E-27 

2.75E-25 

AGw 

Ames 

yes 

Eqn. (5) 
w/Ea=75.2 
kj/mole 
bond 

Girshick 

Girshick 

n/a 

4.09E-07 

ANoX 

Ames 

No 

Krestinin 

Ea=134 

kj/mole 

Krestinin 

Girshick 

n/a 

1.00E-03 


Table 3 Nickel catalyst production of SWNTs as indicated by C0 2 mole fraction for ten combinations of rate coefficients 
and dimer nucleation 


Ni Cluster Model T=1300 K P=30 atm 


Case 

Designation 

Carbonyl 
Rates Ni(CO) x 
x=2,4 

NiCO + Ni n 

NiCO+NiCO 

NiCO 

Dissociation 

Ni 2 Dissociation 

Ni n Eraporation 
8>n>2 

Ni n 

Evaporati 
on n>8 

C0 2 

Produced 
at 0.1 s 
No Ni+Ni 
Nucleatio 
n 

C0 2 

Produced 
at 0.1 s 
Ni+Ni 
Nucleatio 
n 4E15 

1 

AWNiWOldFe 

CO 

Sunderlin 

yes 

0 

Ames Fe 
Ea=24 
kJ/mole 

Ames from 203 
kJ/mole bond 

Girshick from 203 
kJ/mole bond 

Girshick 
from 203 
kJ/mole 
bond 

1.23E-03 

4.80E-04 

2 

AWNiOldFeC 

ONiCONiCO 

Sunderlin 

yes 

Krestinin Fe- 
rate 

Ames Fe 
Ea=24 
kJ/mole 

Ames from 203 
kJ/mole bond 

Girshick from 203 
kJ/mole bond 

Girshick 
from 203 
kJ/mole 
bond 

1.23E-03 

4.80E-04 

3 

AWNiS 

Sunderlin 

yes 

0 

Sunderlin 

Ea=170 

kJ/mole 

Ames from 203 
kJ/mole bond 

Girshick from 203 
kJ/mole bond 

Girshick 
from 203 
kJ/mole 
bond 

n/a 

5.37E-15 

4 

AGNiOldFe2 

Sunderlin 

yes 

0 

Ames Fe 
Ea=24 
kJ/mole 

Girshick from 203 
kJ/mole bond 

Girshick from 203 
kJ/mole bond 

Girshick 
from 203 
kJ/mole 
bond 

1.47E-06 

4.12E-04 

5 

AGNoXNiNoN 

uNoNiCONiC 

0 

Sunderlin 

No 

0 

Ames Fe 
Ea=24 
kJ/mole 

Ames Fe Ea=134 
kJ/mole 

Girshick from 203 
kJ/mole bond 

Girshick 
from 203 
kJ/mole 
bond 

7.72E-25 

n/a 

6 

AWNi 

Sunderlin 

yes 

0 

Sunderlin 

Ea=170 

kJ/mole 

Girshick from 203 
kJ/mole bond 

Girshick from 203 
kJ/mole bond 

Girshick 
from 203 
kJ/mole 
bond 

n/a 

9.92E-25 

7 

ADNiOldFe2 

Sunderlin 

yes 

0 

Ames Fe 
Ea=24 
kJ/mole 

Ames Fe Ea=134 
kJ/mole 

Girshick Fe 

Girshick 

Fe 

3.89E-06 

9.72E-04 

8 

ADNi 

Sunderlin 

yes 

0 

Sunderlin 

Ea=170 

kJ/mole 

Ames Fe Ea=134 
kJ/mole 

Ames Fe 

Girshick 

Fe 

1.19E-19 

1.19E-19 

9 

ADNiCONiCO 

Sunderlin 

yes 

Ames Fe-rate 

Sunderlin 

Ea=170 

kJ/mole 

Ames Fe Ea=134 
kJ/mole 

Ames Fe 

Girshick 

Fe 

3.27E-17 

3.27E-17 

10 

ADNoNuNoNi 

CONiCONoX 

Sunderlin 

no 

0 

Sunderlin 

Ea=170 

kJ/mole 

Ames Fe Ea=134 
kJ/mole 

Ames Fe 

Girshick 

Fe 

3.72E-25 

n/a 


cm /s/mol 


1.E+13 

1.E+12 

1.E+11 

1.E+10 

1.E+09 

1.E+08 

1.E+07 

1.E+06 

1.E+05 

1.E+04 

1.E+03 

1.E+02 

1.E+01 

1.E+00 

1.E-01 
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Fig. 1 Reaction rate coefficients at 1 300 K for iron clusters used in various reaction schemes. Also shown 
the legend is the amount of C0 2 produced after 0. 1 seconds in the various models at 30 atm. 


